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1. Introduction 

There is a long and successful history of calculating the renormalized vacuum 
polarization and stress-energy tensor on spherically-symmetric black holes (see P, El 
El m El El H IB]), however, in the astrophysically significant case of the Kerr- Newman 
black hole, such calculations have proved intractable. This is due to the fact that 
all the calculations in the spherically symmetric case have relied heavily on specific 
properties of the Legendre functions (mainly the Legendre Addition Theorem and the 
Heine Identity) and also the applicability of the Watson-Sommerfeld Formula to obtain 
a mode-sum expression that is numerically tractable. In the axially symmetric case, 
however, addition theorems for the angular functions do not exist nor do we have an 
analog of the Heine identity for the angular functions that arise. We will show, in the 
framework of the Schwarzschild black hole threaded by a cosmic string, that we can 
proceed in the axially symmetric case by obtaining useful summation formulae based 
on the Hadamard structure of the Green function without using specific properties of 
the mode-functions themselves. In this paper we will use this approach to obtain an 
analytic expression for the renormalized vacuum polarization on the horizon of such a 
space-time. 

It should be noted that the result presented here is a generalization of a result 
given by Davies and Sahni [10] where the authors calculate the vacuum polarization for 
a restricted and discrete set of values for the azimuthal deficit. 

2. The Green's Function 

The cosmic string is modeled by introducing an azimuthal deficit parameter, a, into the 
standard Schwarzschild metric. We may describe this in coordinates (t, r, 9, 0), where 
is periodic with period so we may take G [0,27ra), in which the line element is 
given by 

ds^ = -(1 - 2M/r)dt^ + (1 - 2M/r)"Mr2 + rM^^ + sin^ (1) 

For GUT scale cosmic strings a = 1 — 4yU where ^ is the mass per unit length of the 
string so we shall assume < a < 1. We can alternatively define a new azimuthal 
coordinate by 

= a0 (2) 

so that is periodic with period 27r and we may take G [0,27r). In coordinates 
(t, r, 6', 0), the line element is given by 

ds^ = -(1 - 2M/r)dt^ + (1 - 2M/r)-Mr2 + rM^^ ^ ^2^2 ^^^^2 q^^2 _ 

We shall consider a massless, minimally coupled scalar field, (p, in the Hartle- 
Hawking vacuum state. Since this is a thermal state, it is convenient to work with 
the Euclidean Green's function, performing a Wick rotation of the temporal coordinate 
t ^ —IT and eliminating the conical singularity at r = 2M by making r periodic with 
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period 27t/h where k, = 1/(4M) is the surface gravity of the black hole. The massless, 
minimally coupled scalar field, satisfies the homogenous wave-equation 

□¥.(r,r,^,0) = O, (4) 
which can be solved by a separation of variables by writing 

<^{t, r, e, (f)) ~ e*"''"+^'^'^P(^)i?(r) (5) 
where P{0) is regular and satisfies 

while R{r) satisfies 

(^(^2 _ 2Mr)^ - A(A + 1) - I^(^) = 0. (7) 

I dr ^ 'dr ^ > - 2Mr i ^ ^ ^ ^ 

The A(A + 1) term arises as the separation constant. The choice of A is arbitrary for 
if to satisfy the wave equation but requires a specific choice in order for the mode- 
function P{6) to satisfy the boundary conditions of regularity on the poles. In the 
Schwarzschild case (no cosmic string, a = 1), regularity on the poles means that A = /, 
i.e the separation constant is /(/ -|- 1). In the cosmic string case, the appropriate choice 
of A that guarantees regularity of the angular functions on the poles is 

A = / — |m| + \m\/a. (8) 

With this choice of A, the angular function is the Legendre function of both non-integer 
order and non- integer degree, mz.. 

It can be shown that these angular functions satisfy the following normalization 
condition, 

^rJi:'i;vi/.(cos ^)p;ir:/:i^i/ jcos e)d{cos e) = 

2 r(A-|m|/a + l) 
(2A + l)r(A+|m|/« + !)"'■ ^ ' 

The periodicity of the Green's function with respect to (r — r') and {(f) — cf)') with 
periodicity 2-n / k and 27r, respectively, combined with Eq.(lO) allow us to write the 
mode-sum expression for the Green's function as 

^^^'^^ 87r2 2^ ^ 2^ 2^^^^ + ^^r(A-|m|/« + l) 

n=— CO m=-oo l=\m\ 

p-l"^l/"(cose)P,-l"^l/"(cose')XnA(r,rO, (11) 
where Xnx{i^,i^') satisfies the inhomogeneous equation, 

- 2^-)^ - MA + 1) - ;:^^}xnA(r, r') = -^^(r - r'). (12) 
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It is convenient to write the radial equation in terms of a new radial variable 
?7 = r/M — 1, the radial equation then reads 

- '^^) - ^ 1) - i^i^-C'"') - -^^(" - "')■ (^«) 

where we have used the fact that k = 1/4M. For n = 0, the two solutions of the 
homogeneous equation are the Legendre functions of the first and second kind. For 
n 7^ 0, the homogeneous equation cannot be solved in terms of known functions and 
must be solved numerically. We denote the two solutions that are regular on the horizon 
and infinity (or some outer boundary) hy p^^\r)) and q^xKv)^ respectively. A near-horizon 
Frobenius analysis for n ^ shows that the indicial exponent is ±|?t,|/2, and so we have 
the following asymptotic forms: 

pi'Kv) - iv - ip^' v^i, 

Defining the normalizations by these asymptotic forms and using the Wronskian 
conditions one can obtain the appropriate normalization of the Green's function: 

^Px{v<)Qx{v>) n = 0, 

Xnxiv,v') = ^ " 1 I , I , (15) 

-Px{v<)Qx{v>) n=^0. 



2\n\aM 



3. Renormalization 



From the asymptotic forms (14), it is clear that taking one point on the horizon (with 
the other outside) means that all the n ^ modes vanish. The unrenormalized Green's 
function then reduces to 

Gix x') = ^ V e™(^-*') V (2A + n HA + H/g + l) 

^^"^'"^^ 327r2M2« 2^ ^ Z^^^^ + ^^r(A-|m|/« + l) 

m=-oo l=\m\ VII/ / 

P;:l™l/"(cos^^)P^"l™l/°(cos^^')QA(r/). (16) 

In the absence of the cosmic string, i.e. a = 1, this sum can be done using a combination 
of the Legendre Addition Theorem and the Heine Identity, yielding a closed form 
expression [1], that, of course, diverges in the coincidence limit. 

To renormalize we employ the Christensen-DeWitt [9] point-separation method 
subtracting the geometrical singular part of the Green's function for small separations in 
the coordinates. This method relies on the universal Hadamard singularity structure of 
the Green's function which ensures that once we have subtracted the geometric singular 
part, 'f//cr -|- Vlna\ the coincidence limit of the remainder W is finite. To calculate 
{0'^)ren in the massless Ricci-flat case, the only subtraction term needed is 

Gs^na{x,x') = (17) 
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where s is the geodesic distance between x and x' . Having aheady taken one point on 
the horizon we choose to separate radially, placing x' at 2M + e; the geodesic distance 
is then given by 



s(2M,2M + e) 



(2Me)^/2 2 + l/3(e/2M) + 0{t' 



(18) 



,2M (1 - 2M/r' 

Note that this behaves like (Ax)^/^ rather that the usual Ax since the metric is singular 
at r = 2M so g^''"''' = 0(e), but this is simply a coordinate singularity and in no way 
affects the validity of the geometrical subtraction. The subtraction terms up to 0(e) 
are 

1 1 



G 



sing 



327r2Me IQ2ti^NP 
Taking the partial coincidence limit 



+ 0(f 



(19) 



6 —7- 6' in (16), subtracting the 
0) gives us the 



renormalization terms (19) and taking the limit as — t- 1 (e — 
renormalized expression for the vacuum polarization. However, we now face our 
fundamental challenge as the limit cannot be taken in a meaningful way in the forms 
given. In order to take the limit, we must either invent a way to write the divergent 
terms as an appropriate mode-sum and do a mode-by-mode subtraction or alternatively, 
we must attempt to perform the sum in (16) so that we can write the Green's function 
in closed form. In [TT] we follow the former route to obtain numerical values of (<^^)ren 
outside our black hole. However, when possible, the latter is clearly preferred since it 
can give us a simple closed-form solution for the renormalized vacuum polarization. In 
the next section, we derive a formula that will allow us to perform the sum (16). 



4. Generalized Heine Identity 

We recall that for real arguments, the Heine Identity for the Legendre Functions is [T2] 

oo 

^(2/ + mmQiio = j^-^y (20) 

valid for |\E'| < \C^\. Writing \E' = cos 7 = cos 9 cos 9' + sin ^ sin ^' cos (0 — (p'), and using 
the Legendre Addition Theorem, 

PKCOS7) = ^ e^-^'^-^'^fj^Pricos9)Pr{cos9'), (21) 
^ — ^ (/ + ml! 

m=—l 

we can re-write the Heine Identity as 

E E e^'-^'^-^'Hs/ + i)|^^pr(cos^^)pr(cos^^')QKC) = 7-^ — (22) 

Uii^-i (^ + ^)' (C-COS7) 

Below we shall generalize this form of the Heine Identity to the non-integer Legendre 
functions that arise in the mode-sum of a cosmic string space-time. It is this form 
rather than a generalization of (20) that is useful to us since the fundamentaly axially- 
symmetric nature of the Green's function for a 7^ 1 means that its radial part will 
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always depend on m (through A), and therefore the m-sum over the angular functions 
can never be done independently of the radial part. 

In the Appendix, we have derived other identities involving Legendre Functions of 
non-integer order and/or degree. We have done so since, although they are not directly 
useful here, they turn out to be useful in other contexts, for example, determining the 
Hadamard singularity structure of a 3D Green's function on a cosmic string black hole 
space-time [TTj . 

Returning to the problem at hand, we want to construct a generalization of (22). 
In addition, we would prefer to do this in such a way that the derivation does not rely 
on specific properties of the Legendre functions themselves, since ultimately we wish 
to apply the same techniques in Kerr space-time. There is a very natural way to do 
this by equating equivalent expressions for the same Green's function on a conveniently 
chosen space-time. In our case, the most appropriate case to consider is that of an 
infinite cosmic string in otherwise fiat space-time. Of course, normally in this case one 
would exploit translational invariance along the string and work in cylindrical polar 
coordinates but for our purposes we choose to work in spherical polar coordinates, so 

ds^ = dr^ + dr^ + rM^^ + aV^ sin^ ^d^^ (23) 

where (f) is periodic with period 27r. The Green's function for a minimally coupled scalar 
field on this space-time has a closed-form solution, which can be obtained, for example, 
by separating in cylindrical polar coordinates [13] 



G(x, x') = — ^^Mx/cy) ^ .24) 

Svr^a pp' sinh x(cosh(x/ a) — cos(0 — 0')) ' 

where 

(r — r')^ + + r'^ — 2rr'cos6'cos6'' , , 

coshv = (25) 

^ 2pp' ^ ^ 

and p = rsin^^ is the cylindrical polar coordinate. 

One can also obtain the spherical polar mode-sum expression for this Green's 

function: 

•^^ m=-oo l=\m\ ^ ' ' 

p-l'"l/"(cose)P,-l'"l/"(cose')^^A(r,r'), (26) 
where Rujx{r,r') is the ID Green's function satisfying 

The solution of this equation can be given in terms of Modified Bessel Functions, 

(28) 

where / and K are the Modified Bessel Functions of the first and second kinds, 
respectively. The integral over u can now be performed using the identity [1] 



POO ^ ^ 

/ cosa;(r - r )/A+i/2(c^r<)irA+i/2(c^r>)du; = ;ni/2 ^^(0 (29) 
Jo ) 
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valid for A > — 1, where Qx is the Legendre function of the second kind and 

C = . (30) 

Thus, the spherical polar mode-sum expression for the Green's Function may be written 
as 

r(A + |m|/a + l) 



G(x,x') = — ^ y e'-(<^-*') y (2A + 1); 



■^n'^arr' ^ ^ FfA— |m|/a + l) 

m=-oo l=\m\ VII/ ) 

P;I'"I/"(cos^)P;I'"I/"(cos0')Qa(C)- (31) 

We now equate this Green's Function with the equivalent closed form expression 
54|) to obtain the following generalized Heine Identity: 



sinh(x/Q;) 



F(A- |m|/« + l) 

m=-oo l=\m\ \ I II I 



(32) 



sin6'sin6''sinhx(cosh(x/a) — cos(0 — 0')) 
where 

coshx=^^ 7, 7 — ^ C = COS7 + (coshx — cos(0 — 0')) sin6' sin6''. (33) 

sm d sm d 

This identity, valid for C > | cos7|, is completely analagous to (22) and reduces to it in 
the a — )■ 1 limit. 



5. Calculation of (i^^). 



Returning now to our calculation of the vacuum polarization on the horizon of the 
Schwarzschild black hole threaded by a cosmic string, we see that our generalized Heine 
Identity (32) has precisely the right form to enable us to perform the mode-sum in the 
Green's function expression ( 11 ). The Green's function, where one point has been taken 

sinh(x/«) 



to lie on the horizon, may now be written as 

1 



' 327r'^M'^a sin 6 sin 6' sinh x(cosh(x/ a) — cos(0 — 0')) 
For purely radial separation, we take 0' — 0, ^' — 6*, to yield 

1 sinh(x/a) 



G(x,^;0 

where we now have 

coshx = 



327r^M2a sin^ Q sinhx(cosh(x/a;) 
?7 — cos^ Q 



sin^^ 



1 + 



Msin^^' 



(34) 



(35) 



(36) 



on taking r' = 2M + e <(=^ r] = e/M+1. From Eqs.(19) and (35), the renormalized 
vacuum polarization is given by 



2 \ horizon 
ren 



lim 

€-S>0 



lim 

e-S>0 



G(X,0) - Gsingix,0) 

sinh(x/a 



327r2M2asin2^sinhx(cosh(x/a) - 1) 32n^Me 192tt^M^ 



(37) 
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It is now convenient first expand in terms of x 

sinli(x/a) 2a / 1 a 



tlien 



)+0{x'). (38) 



so we obtain 



sinli x(cosli(x/a) — 1) X^ 3 

e 1 Msin^^ ^ ^, i, 

sinli(x/a) aMsin^^ 1 2 



1 / w / X X + — (1 - a') + 0(e^). (40) 

sinfix(cosli(x/«) - 1) e 6a ^ ) \ ) \ ) 

Substituting ( [40| ) into ( [37| ) and taking tfie limit, we arrive at tlie renormalized vacuum 
polarization on the horizon, 

<^^>""'°"-l92^(l + S)' 

The first term here is simply the Candelas result [T], which, of course we recover 
when a — )■ 1, and the second term represents the contribution due to the presence of the 
cosmic string. The result (41) has in fact been calculated by Davies and Sahni [10] but 
for a restricted set of a values such that (1/a) is an integer. For this restricted class of 
values, the cosmic string solution is obtained as a sum of (1/a) images of the solution in 
the absence of a string. The method of images cannot be applied for general azimuthal 
deficits, however, so our method presents a generalization and independent derivation 
of the result of Davies and Sahni. 

We have plotted the vacuum polarization in Figure [T] for a range of values of a. We 
see that the presence of the string increases ('^^)re"^°" everywhere. In the equatorial 
plane it is increased by a factor of 1/a^ from its value in the absence of a string. As we 
approach the axis (<^^)re"^°" diverges in a non-integrable manner as 

-.^0. (42) 

as one would expect a distance (2Msin^) from a flat-space cosmic string. The region 
over which the string dominates may be characterised by the range of cos 9 (near = 0) 
at which (<^^)re"^°"' exceeds twice its value in the equatorial plane, this is given by 

1 - (cos^^)2 = 1 - ^2)1/2 = (!-«) + 0((1 - af)). (43) 



6. Conclusions 



We have considered a massless, minimally coupled quantum scalar field on a 
Schwarzschild black hole threaded by an infinite straight cosmic string. We have 
obtained a simple analytic expression for the vacuum polarization on the black hole 
horizon for this field in the Hartle-Hawking vacuum state, which is a generalization of the 
Candelas \[\ result. In order to do this, we have derived a very useful summation formula 
involving non-integer Legendre functions, which we have called a generalized Heine 
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Figure 1. {<p'^)ren^°" (ill units M = 1) as a function of cos 6 for a range of values of a. 
For a ^ I, {'p'^)ren^°'^ diverges at the poles since there is a curvature singularity there 
due to the presence of the string. {(p^)rtn^°^ increases everywhere as a decreases. 



Identity given its analogy with the standard Heine Identity involving Legendre functions 
of integer degree. The derivation of this formula involved equating equivalent expressions 
for the same Green's function on a conveniently chosen space-time. Furthermore, the 
derivation did not rely on specific properties of the Legendre functions. This idea can be 
used to derive many useful summation formulae and addition theorems for non-integer 
Legendre functions, some of which are presented in the Appendix. Elsewhere [TT], we 
have used this idea to obtain an expression for the vacuum polarization on the entire 
exterior region of the space-time. We believe this framework will prove useful in the 
other axially symmetric black hole calculations, most importantly, the Kerr-Newman 
black hole. 



Appendix 

In this Appendix, we shall derive some interesting summation formulae for non-integer 
Legendre functions based on considering equivalent forms of the 3D Green's function on 
flat space threaded by a cosmic string 

ds^ = dr^ + r^d^^ + aV^ sin^ ^d^l (44) 
The 3D Green's function for Laplace's equation satisfying 

V2G3(x, x') = -^?3(x)-^/2<5(x - x') (45) 
can be written in spherical polar coordinates as the mode sum: 



Ana ^ ^ r(A - |m|/a + 1) 

m,=-oo l=\m\ ^ ' " ' 



p-l"^l/"(cose)P~l'"l/"(cos^^')(^<)\^>)"^^'''^- (46) 
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Note that this corresponds to the n = term of the 4D Green function in the 
corresponding ultrastatic {gu = ^1) space-time, but that the corresponding statement is 
not true in the merely static Schwarzschild cosmic string space-time. We can also write 
this same Green's function in several other forms and in several coordinate systems, for 
example, performing the mode decomposition in cylindrical polar coordinates, p = r sin 6 
and z = rcos9, we have 

1 °^ poo 

Gs{^,^') = — J2 e*™^'^"'^'^ / cosk{z~z')I\m\/a{kp^)Kim\/a{kpy)dk. (47) 
The integral over k may further be performed to yield 



^ ' m=— OO ^ 

Another form of this Green's function is obtained from General Axi-Symmetric Potential 
Theory [l5l[T6]: 



47r^a 

m=— oo 



(pp')|m|/Q sin2|m|/a ^ 



/o [(2-20' + P' + P"-2pp'cos*](H/-+i/2) • ^ > 

These expressions are valid for all a G IR+; equating any one of these expressions with 
another gives us a formula for the functions involved. Some of these results are well 
known and can be found in volumes such as Gradsteyn & Rhyzik [ll], nevertheless, we 
have non-obvious, important results, for example, equating (46) with (48) and taking 
r — )■ r' we obtain 



oo 



V r(A + \m\/a + Ij p-H/^/^^gmp-H/g/^Qgg/x 
^ r(A-|m|/a + l)^^ icosfyjf^, (cos^j 

l=\m\ ^ II/'; 



1 ^ /1-cos^cos^^, 

T:[smU smU')^i'^ \ smb'smfc' / 



There are also a set of formulae that hold for a more restrictive set of a values. 
Linet [TTj has shown that, when a > 1/2, the Green's function may be written as 

G3(x,x') = — + — / -F„(n,0-0OdH (51) 

where 

a = [p2 + p'2 _ 2pp'cosa((/) - </>')]'/' 

= [p2 + p'2 ^ _ ^ 2pp' cosh m] (52) 

and 

P'i:'^' *) = f "('^ - '^/") + , ^ . (53) 

cosh(-u/aj — cos(W — vr/aj cosh('u/aj — cos(\l/ + vr/aj 
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Equating this to ([46]) for a > 1/2 and taking r — )■ r', we find 



a r(A — |m|/a + 1 

m=— cxD i=|m| 



1 ^ 1 1 /"^ _FUw,0-£^ 

+ ^7TT— TT— TiTNTT^ / r.„.t, ^ , .^.u^.n/^ ^^ (54) 



[2(l-cosa7)]V2 27ra (2sin^sin^')^/^ Jo [cosli^ + cosh u] 1/2 
where 

cos 07 = COS 9 cos 6*' + sin 9 sin 6*' cos a(0 — 0') 

, ^ 1 — cos ^ cos ^' 

= ■ a ■ Of ■ (^^) 

sm 9 sm 6' 

This formula has the benefit of relating a 3D mode-sum expression to an expression that 

isolates the 3D Hadamard singularity structure (the first term on the RHS) together 

with a regular term that depends on the boundary conditions. We have used this result 

to prove the regularity of the mode-sum that arises in the calculation of the vacuum 

polarization on the exterior region of the spacetime [H]. 

Another form of the Green's function can be given in terms of the half-integer 

Legendre functions by considering the Green's function in toroidal coordinates, which 

are related to the cartesian coordinates by 

sinh /i cos (f) sinh fi sin sin rj 

X = — y = — z = — , (56) 

cosh /i — cos rj cosh // — cos rj cosh /i — cos rj 

where the ranges of the coordinates are O</i<oo,O<r7<27r,O<0<27r. The 

mode-sum form of the Green's function obtained by separating in these coordinates is 

00 00 

G3(x,x') = -^ y e*"^(<^-*') y e*"('^-'''H(cosh;U- COS r/) (cosh COS V)]^/^ 

m=— 00 71=— 00 

Equating this to (48) gives us the following addition theorem 
E e-'-'-'')e""l/.. J(" + H/° + ;/^) p,;lH/°(coBh^,)(3.:H^(cosh^,) 



n=— 00 



r(n - \m\la + 1/2)' 

Q|m|/Q-1/2(X) 



[sinh /i sinh /i' (cosh /i — cos r^) (cosh /i' — cosr]')]^/^ 
The final coordinate system we shall investigate is prolate spheroidal coordinates 
which may trivially extended to the oblate spheroidal case. The relationship between 
prolate spheroidal coordinates and cartesian coordinates is given by 

X = sinh cr sin 6' cos y = sinh a sin ^ sin 2; = cosher cos ^. (59) 

where the ranges of the coordinates are O<cr<oo,O<^^<7r,O<0<27r. The mode 
form of the Green's function obtained by separating in these coordinates is 

nj^ - — V P^™('^-'^') V r9A I 1 vlm|/«^i r(A+ |m|/a + 1)^ 
^3lx,xj-^^^ e 2^i2A + ije r(A-|m|/a + l)2 

m=-oo l=\m\ \ I 1/ I 

p-l'"l/"(cos0)P^"l"l/"(cos0')^A"'™'^"(cosha<)g-l'"l/°(cosha>) (60) 



Schwarzs child Black Hole Threaded by a Cosmic String 



12 



We shall write down only one possible summation formula here, by equating the Green's 
function expression above to (48), we obtain a summation formula for a product of four 
Legendre functions: 

i=\m\ \m\ia + i) 

p-H/"(cosha<)Q-H/.(,^3ha>) = - ^=#^^^P^^== (61) 

Tra V smh a smh a' sm 6 sm U' 

where 

cosh^ a + cosh^ a' — sir? 9 — sii? 9' — 2 cosh a cosh a' cos 9 cos ^' . ^, 

X ~ (62) 

2 sinh a sinh a' sin 9 sin 6'' 

There are many more summation and addition formulae that can be derived in this way. 

Those presented here were either of particular use to us elsewhere or are to the best of 

our knowledge previously unpublished. 
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